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SLAMCODE Computer Program 


The program SLAMCODE developed in this study and 
described in this Volume is a large Finite Element Method 
(FEM) code capable of dynamic analysis of transient problems. 
It was adapted from the original SLAMCODE written by C. aks 
Costantino (Costantino, 1968) by incorporating the implicit 
Newmark B-integration (Newmark, 1959) and by constraining 
the free surface to be stress free. The present code also 
differs from the original by not permitting plastic elements. 
It calculates solutions for axisymmetric, plane stress or 
plane strain problems, as in the original code. It is designed 
to run on the IBM 360 family of computers under 0S/360, Being 
written in Fortran, it can be easily adapted to other com- 
puter systems. In fact, it was written with this requirement in 
mind by minimizing the number of system-dependent Fortran 
peculiarities. 

The Fortran source cards are listed in Appendix A, The 
Link Editor control cards are listed in Appendix B. These 
control cards are appended to the object module for execution 
of the program. Appendix C contains examples of Job Control 
Language cards for the Fortran H computer using either two or 
four 9-track tape drives. Tables 1-4 describe the operation of 
the program and Figure 1 is a block diagram of the overlay 
structure. Finally, detailed descriptions of the'input para- 
‘meters and operation of the program can be gleaned by reading 


the comments in the source listing. All input data is described 
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where it is read in the program listing. 

The bandwidths of the diagonal matrices are cortrolled 
by the numbering of the nodal points. The bandwidth minimiza- 
tion algorithms listed in Appendix A are those used by 
Costantino, (PATH) and (MINI). A third algorithn written by 
Rosen (1968) and adapted for SLAMC@DE is listed in Appendix Al, 
(MINBND) . The listing of the original Rosen algorithm con- 
tained some typographical errors which were corrected for 
(MINBND), but the comment cards which we wave not reproduced 
in Appendix Al are very useful. (PATH) is much faster than 
(MINBND) but its effectiveness depends upon the selection of 
starting node points. (MINI) has not decreased the bandwidth 
for any of our element arrays. 

In general, large problems are solved by several runs of 
the program, all output being saved on magnetic tape or 
cataloged disk files. A wide range in input/output (1/9) 
flexibility is available with 9S/360 and this should be exploited 
to maximize program efficiency. For example, tape units should 
only be used in preference to disks when the time reduction made 
possible by saving intermediate results exceeds the time lost by 
the tape drives. 

The program in its present form is I/P bound, i.e., a large 
amount of the computer time is spent spinning tapes and disks. 
In a typical 360/67 run there are approximately the Same number 


of I/f seconds as there are CPU seconds. This performance 


might be improved on larger machines where I/Q@ buffers and 
program dimensions could be increased. The 360/67 version 
runs in a 280K byte region of core memory. 

The number of I/® devices required depends on the 
operating system. The 360/67 system can assign many files to 
one device. A general rule of thumb for a system with 280K 
bytes of available core memory is that the program requires at 
least four tape drives and a hundred cylinders of disk file 


space. 
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